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ABSTRACT 

In  this  analysis  a  generalized  method  is  developed 
for  determining  the  optimum  stage  weight  distribution 
for  multistage  rockets.  Inclusion  of  the  variations  in 
structural  factors  with  stage  weights  in  the  optimization 
process  is  shown  to  lead  to  a  more  generalized  set  of 
optimum  conditions.  Expression  of  all  rocket  weight 
parameters  in  terms  of  the  stage  weights  allows  for 
convenient  optimization  as  well  as  for  a  comparison  with 
previous  optimization  methods. 

This  approach  permits  improved  optimum  design 
over  existing  methods  for  maximizing  payload  ratios 
for  given  ranges  and  for  maximizing  ranges  for  given 
payioad  ratios.  An  evaluation  of  previous  methods  is 
included  for  comparison  purposes,  and  the  limitations 
of  these  previous  methods  are  discussed. 


nt 


STL./  TN- 60 -0000 -090  36 
Page  iii 


CONTENTS 

Page 

ABSTRACT .  ii 

I.  INTRODUCTION .  1 

NOMENCLATURE  .  . .  2 

II.  OPTIMIZATION  PARAMETERS  .  3 

A.  Performance  Parameter* .  3 

B.  Structural  Factor  Parameter* .  4 

C.  Use  of  Stage-Weight  Parameter* .  5 

III.  THE  GENERAL  OPTIMIZATION .  7 

A.  Lagrangian  Multiplier  Technique .  7 

B.  Application  of  the  Lagrangian  Multiplier  Technique  to 

Optimum  Rocket  Design .  8 

C.  Criticism  of  Previous  Method* .  8 

D.  Reduction  of  the  Optimization  Equations  to  Simpler  Form  .  .  11 

IV.  APPLICATION  TO  OPTIMUM  THREE-STAGE  ROCKET  DESIGN  .  12 

V.  APPLICATION  TO  OPTIMIZING  TWO-STAGE  ROCKET  DESIGN  .  15 

VI.  APPLICATION  TO  N-STAGE  ROCKET  OPTIMIZATION .  16 

REFERENCES .  17 

APPENDIX .  18 


STL/  TN -60 -0000 -09036 

Page  1 


I.  INTRODUCTION 

In  recent  years,  the  optimization  of  multistaged  rockets  has  received 
considerable  attention.  Since  the  performance  of  missiles  and  space  vehrcles 
is  sensitive  to  small  changes  in  design,  optimization  procedures  are  of  great 
importance.  In  References  1  through  11  methods  were  developed  for  deter¬ 
mining  "optimum"  stage  mass  ratio  distributions.  None  of  these  methods, 
however,  allowed  for  variations  of  structural  factor  with  stage  weight.  Conse¬ 
quently,  the  referenced  methods  do  not  yield  truly  optimum  stage  mass  ratio 
distributions. 

The  purpose  of  this  paper  is  first,  to  derive,  in  terms  of  the  stage  weights, 
a  more  general  design  optimization  that  allows  for  variations  of  structural 
factors  with  stage  weights,  and  second,  to  evaluate  the  limitations  of  previous 
design  criteria. 
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NOMENCLATURE 


r . 
J 

R 


D  and  B 


Payload  ratio,  as  defined  by  Equation  (1) 

Payload  weight 

Grot*  vehicle  weight 

Burnout  velocity 

Specific  impulse  of  the  j**1  sta^e 

Acceleration  of  gravity 

Mass  ratio  of  j**1  stage 

Range 

Empirical  parameters  for  range  Equation  (3) 

Lumped  velocity  requirement  term 

Structural  factor  of  i^  stage 

Weight  of  the  i^  stage 

Propellant  weight  of  the  ith  stage 

Gross  weight  of  ith  stage  as  defined  by  equation  (11) 

Total  number  of  stages 

Partial  derivative  of  f  with  respect  to  W.,  keeping  <r.  fixed 
Partial  derivative  of  f  with  respect  to  <r.,  keeping  W  fixed 
Burnout  weight  of  the  ith  stage,  as  defined  by  Equations  (12)  and  (13) 
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II.  OPTIMIZATION  PARAMETERS 


A.  Performance  Parameters 

The  performance  capability  of  a  multistage  rocket  vehicle  car.  be  described 
by  two  equation*. 

wd> 

p  ■  -wr  (i) 

la 


N 

=  E  I.  g  In  r. 
bo  j=i  J  J 


-  6  V 


(2) 


where  6  V  represents  the  velocity  losses  associated  with  gravity  and  drag.  The 

drag  losses  are  primarily  dependent  upon  the  initial  thrust-to-weight  ratio,  N  . 
and  on  the  quantity,  VC^A.  ° 

Equation  (2)  can  be  rewritten  in  terms  of  range,  R,  for  a  ballistic  missile: 
(Reference  12) 


R 


N 

i n 


WB 


-  1 


(3) 


where  B  is  very  insensitive  to  changes  in  Nq  and  W^/CpA,  while  the 
parameter  D  is  fairly  sensitive  to  such  change*. 

Let 


and  the  theoretical  velocity, 


(4) 


♦  6  V 


(5) 
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Then  from  Eqs.  (2),  (3),  and  (4), 


or 


(6) 


(7) 


Thus,  for  a  given  initial  thrust-to-weight  ratio  and  mission,  the  velocity 
requirements  may  be  lumped  into  a  fixed  single  term,  The  range  equation 
already  provides  for  velocity  losses  in  the  empirical  constants  D  and  B  used 
for  any  particular  configuration. 

When  there  are  at  least  two  stages  and  when  all  the  specific  impulses  are 
known,  P  and  f  do  not  uniquely  define  the  mass  ratios  of  the  various  stages. 
Consequently,  proper  selection  (optimization)  of  mass  ratios  for  either  maxi¬ 
mum  payload  at  a  given  range  or  maximum  range  at  a  given  payload  ratio  is 
required. 

B.  Structural  Factor  Parameters 

The  structural  factor,  cr,  for  the  stage  is  given  by: 


<r. 

l 


wi- 


W 


(8) 


Expressed  in  terms  of  the  weight  of  the  i**1  stage,  the  following  scaling  laws 
are  assumed  to  hold: 


9 


i 


n  -1 

ciw, 


(9) 


where  C.  and  m  are  empirical  constants  for  each  stage  subject  to  the  selection 
of  propellant  feed  systems,  auxiliary  systems,  etc. 


STL/TN-60-0000-09C36 

Page  5 


C,  Use  of  Stage-Weight  Parameter* 

The  mass  ratio  of  the  i**1  stage  can  be  defined  by: 


r.  =  - p-r 

i  ^TTT, 

bo 


W'*'  -  W  . 
o  pi 


(10) 


where 


(i)  N 

Wu'  =  Z.  W.  +  Wt 
o  )  L 


(11) 


wb»i)  *  *:»  -  %  <u» 

From  Eqs.  (8),  (11),  and  (12),  the  following  useful  relatione  can  be  derived  for 
an  N- stage  rocket: 

W,  =  a.  W,  +  W*2)  =  <r  +  (1  -  OW*2* 

bo  11  o  lo  1  o 


W,  *2)  =  <r,  W,  +  \v  =  <r,  W*2)  +  (l  -  «r_)  W(3) 
bo  2  2  o  2  o  2  o 


WboN>  *  'H  WN  +  WL  -  'N  WiN)  +  "  -  «N»  WL 


Combining  Eqs.  (9J,  (10),  (11)  and  (13), 


N 

S  W.  +  W, 


(13) 


Ci  Wl  ♦  J-fcl  Wi  +  WL 
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Substituting  (11)  into  (1), 


P 


V  WL 

WL 


(15) 


wliere  the  payload  ratio  is  now  expressed  in  terms  of  the  stage  weights.  By 

substituting  (14)  into  (6)  or  17),  f  can  also  be  expressed  in  terms  of  the  stage 
weights: 


♦ 


r 


vs 


i 


N 

Sc  Wi 


L 


+  w, 


w. 


N 


j«kl 


W. 

J 


+  w. 


(16) 


% 
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III.  THE  GENERAL  OPTIMIZATION 

A.  Lagrangian  Multiplier  Technique 

A  necessary  condition  that  a  function  f{Xj,  x^,  ...»  x^)  of  N  variable* 

Xj,  x^ . Xj^  have  a  stationary  value  i*  that 

df  =  Iq-  dxl  +  W~  dx2  *  •  ”  ^  =  0  (17) 

for  all  permissible  values  of  the  differentials  dXj,  dx^, . . . ,  dx^.  If,  however, 
the  N  variables  are  not  independent,  but  are  related  by  another  condition  of  the 
form  ^(Xj,  x^, ....  Xj^)  =  0,  then  the  procedure  of  introducing  the  so-called 
Lagrange  multiplier  may  be  conveniently  employed: 

Now  if  X  is  so  chosen  that 


0f  +  \  8*  .  o 

+  X  J—  -  0 


w  +  . 
sir:  +  x 


(18) 


df 


+  X 


* 


0 


4^  (xi  •  x2 . =  0 

then  the  necessary  condition  for  an  extremum  of  f(Xj,  x^,  ....  x^)  will  be 
satisfied.  The  quantity  X  is  known  as  a  Lagrangian  multiplier. 
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B. 

The  Lagrangian  multiplier  method  may  be  used  to  optimize  P  subject  to  a 
fixed  $  or  vice  versa.  In  fact  the  partial  differential  equations  in  (18)  would 
be  the  same  for  either  optimization;  only  the  constraining  equation  would  be 
different. 

In  general,  the  conditions  (18)  applied  to  optimum  stage  weights  become: 


>plication  of  the  Lagrangian  Multiplier  Technique  to 
itimum  Rocket  Design 


8  P  %  8* 


0 


and  either  $  or  P  =  constant 
C.  Criticism  of  Previous  Methods 

The  optimization  conditions  expressed  by  (19)  guarantee  a  minimum  P  for 
constant  £  or  vice  versa,  since  variation  of  structural  factor  with  stage  weight 
is  included  in  the  optimization  process.  In  previous  methods  (References  1  to 
11),  this  variation  was  not  included.  To  evaluate  these  methods,  0  can  be 
re-expressed  as: 
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Then: 


Equations  (19)  then  become 

f 

I 

i 


I 

t 

(21) 


(22) 


and  either  P  or  $  constant. 

When  no  variation  of  structural  factor  with  stage  weights  is  considered,  the 
corresponding  expressions  for  "optimization”  can  easily  be  shown*  to  be 


* 

See  Appendix 
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OP 

W 


l 


8  P 
WT 


0 


0 


(23, 


8  P 


0 


and  either  P  or  ^  =  constant. 

Since  relations  (22)  are  the  true  optimization  conditions  then,  in  order 
for  relations  (23)  to  be  optimum,  (22)  and  (23)  must  be  compatible.  It  is  clear 
that,  in  order  for  (22)  and  (23)  to  be  compatible,  the  relations 


JLt 

8«r  i 


lw. 


drl 

^7 


0 


(24) 


-  <rN 


0 


must  hold. 


STL/ TN -60-0000-09036 
Page  11 


Since  (24)  cannot  be  valid,  except  in  the  nonrealistic  case  where  the 
structural  factor  does  not  vary  with  stage  weight,  Eqs.  (22)  and  (23)  are  not 
compatible.  Hence,  Eqs.  (23)  do  not  represent  a  true  optimization  criterion 
for  realistic  rocket  design.  The  actual  discrepancy  in  the  design  criterion  in 
using  (23)  rather  than  (22)  will  depend  upon  the  relative  magnitudes  of  the 
terms  in  (22)  that  have  the  coefficient  X.  If  the  second  term  of  coefficient  X 
in  (22)  is  negligible  compared  to  the  first  term  of  coefficient  X,  then  (23) 
represents  a  realistic  optimization.  In  general,  however,  this  is  not  the 
case. 

D.  Reduction  of  the  Optimization  Equations  to  Simpler  Form 
From  Eq.  (15),  for  fixed  W^: 

8  P  3P  8P  _  1  ,, 

JW^  "  aw2  ‘  •••  ~  8W^  ~  U 

Substituting  (25)  into  (19)  and  eliminating  common  terms,  the  optimization 
equations  reduce  to: 


il  _  8  *  _  d± 

awj  "  ~  “  fw 


(26) 


and  either  0  or  P  =  constant 
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IV.  APPLICATION  TO  OPTIMUM  THREE-STAGE  ROCKET  DESIGN 


For  a  three-stage  rocket,  Eqa.  (15)  and  (16)  become 


P  = 


W1  +  W2*  W3  +  WL 

- w: - 


(27) 


and 


(  W.+ 


1 


w,  +  w,  +  w, 


\X1  f 


n 


\^1  "1 
Accordingly, 

<TW^“ 


C,W,‘  +  W,  +  W,  ♦  V/, 


w,  +  w„  +  w. 


n. 


VC2W2  +W3  +  WL 


W,  +  W,  \ 


\I. 


n. 


C3W3  tWV 


i)[<‘  -"i)ciwi‘  +  (»  "M  C1W?  ')  (wz +  w3 +  WL>] 
(W1  +  W2  +  W3  +  W  J  (C1  Wl'  +  W2  +  W3  *  Wl) 


(28) 


(29) 


.  _ »i,(c,w"1''-i)wi _ 

2  (W1+W2  +  W3  +  WJ(C1W?  +  W2<W3  +  Wl) 


* 


f 

I,!(l  -n 


n 


2  (  n2  "M 

*  i  ^  nf  ■»  !  /  tir  ^  ^  | 

L  . 


2)C2W2‘ni.n2C2W2‘  j,  W3 


(W2  +  W3  +  WJ  (C2W22+W3  +  Wl) 


(30) 
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and 


_ _ 

(W1  +  W2  +  Wj  *  wj  (C,  w"1  +  W2  ♦  W3  ♦  wj 

(w2  +  w3  +  wj  (c2  w22  ♦  w3  ♦  wj 

+  »<3[(l  -n3c3W3rl)wL] 

(W3  +  WJ(C3W33  +  Wl) 


(31) 


Substituting  (29),  (30)t  and  (31)  into  the  optimization  equations. 


a  <p  _  a »  _  a  0 

awj  "  WJ  "  aw3  • 


and  factoring  out  0  yields  two  independent  nonlinear  equations: 


!,(■  -n.C.w"1'1) 

(*  -C1W?  ') 


n,  - 1 


1  -  ci  wi 


f  w> 


♦  w„  ♦  w„  ♦  w 


L 


n 


=  I. 


1  "  n2C2W2 


vc1wr  +  w2.w3+wL 

n,-l f  W2  +  W3  +  WL\ 


n. 


\?zw2*w3*wl 


yj 


(32) 
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i,  V- 


new  * 

2  2  2 


'2-^ 

=n 


*-«*** 


*  i. 


1  -  »3c3w, 


n3-li 


V'lJVV^t 

n2 

C2W2  4  ♦  w3  ♦  w. 


W3*WL 


c,w,  5  +  w 


3  3 


(33) 


Equations  (32)  and  (33)  together  with  the  constraint  equation,  f  or  P  = 
constant,  constitute  three  equations  in  terms  of  the  three  unknowns 
W j ,  W2,  and  W^.  For  specified  sets  of  values  of  1^,  I2>  Ij,  Cj,  C2,  Cj, 
nj,  n2,  n?,  and  f  or  P  (depending  upon  which  type  of  optimisation  is 
desired),  the  three  equations  can  be  solved  by  iteration  techniques  to  yield 
optimum  values  of  W^,  and  W2,  and  W^.  By  employing  Equations  (1 1)  and 
(13),  the  other  rocket  parameters  can  also  be  determined. 


\ 
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V.  APPLICATION  TO  OPTIMIZING  TWO-STAGE  ROCKET  DESIGN 

Analogous  to  the  three  stage  optimization,  the  optimization  conditions  for  a 
two  stage  rocket  become: 

at 

*  TTv^ 

(34) 

and  4*  or  P  *  constant 


which  yield: 


(l  -  n,  C,Vl  -  ■) 
‘  (‘  *  ClWin*  '  ‘) 


1  '  C1W1  * 


wi  *  we  ♦  wl  A 

Fl"l“l+  W2  +  WJ 


c2  W.,n2 


WZ  +  WL 

Fz  w2ni  ♦  wL 


(35) 


Equation  (35)  together  with  the  coniraint,  t  or  P  =  constant,  can  be  solved 
by  iteration  to  yield  optimal  values  of  Wj  and  W^. 
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VI.  APPLICATION  TO  N-STAGE  ROCKET  OPTIMIZATION 

If  iterative  techniques  can  be  expanded  to  include  solving  N  simultaneous 
nonlinear  equations,  then  the  process  described  above  can  be  extended  to  any 
number  of  stages.  In  the  case  of  a  large  number  of  stages,  computerized 
random  search  techniques  might  be  employed  to  solve  for  optimum  values  of 
the  stage  weights.  Once  the  optimum  stage  weights  are  obtained,  these  values 
can  be  substituted  into  Eqs.  (11)  and  (13)  to  solve  for  the  other  rocket 
parameters. 
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APPENDIX 


Further  comptrlion  with  prtvitui  methods  writing  4  lor  a  throe  stage 
rocket  in  terms  of  r  ^  and  fg. 


I.  I, 

Wj  t  W2  ♦  Wg  +  WL  \  1  f  W2  +  Wg  +  WL  N  & 


Then 


r^rprwp  Wj  +  wL 


wJ  +  wL 


3  3 


(H) 


8,  |  !,♦<'  -  'i>(WJ  ♦  W3  +  WL) 

'  T^'TWg  <■'  w3  v  ffjjw,  wL) 

88  I  ’  ♦*'!■  *  1,W. 

riq],  *  TSjVfrw^— w3"+  w£j<w3  +  w2  +  w3  ♦  wLj 

I28(l  -  »2)(W}  +  wL) 

*  {sprpmyr  wL)(w2  +  w3  +  wL) 


(W> 


(38) 


and 


Ijtf'j  .  l)Wj 


,  *  l  rw2  ♦  vrp-wLHW7 v  w2  +  w3  +  wj 


*2  '  *>w2 


4  r«r2v?2  ;  w3  vw’-y-(W”T  ’w3  rwj 

13^1  '  VWL 

4  ♦  wL)(w3  ♦  wj 


09) 


Sub*  tut  mg  (37),  (38)  and  (39)  into  the  conditions  from  (23): 


*  &i2  *  ^r|c3 
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result  in: 


1  . 


Wj  +  w3  ♦>  wL 

M'rjW,  +w2  + WJ  +  WL 


=  I. 


1  - 


1  - 


♦  W3  4  WL 
'jI^WjTTjTW' 


w3»wL 
'3  iiT^Tnr-i 


(40) 


Equation  (40)  represents  the  relations  corresponding  to  previous  methods  (1-11) 
for  design  criteria.  By  using  Eqs.  (11)  and  (13),  then  Eqs.  (40),  (32),  and  (33) 
can  be  written  for  comparison  purposes  in  terms  of  the  mass  ratios.  Equation 
(40)  then  becomes 


IjO  -e^)  =  I2U  ‘<r2r2)  c  !3(1  *  r3r3>  (41) 

which  is  the  more  familiar  form  usually  seen  in  the  references  listed  above. 
Equations  (32)  and  (33)  become: 


(1  - 

(1  -  <r,r,)  =  I2(l  -  n2»2r2) 

(42) 

(1  -  frj) 

(1  -  n,  «r  ) 

as 

cr^r 

(1  -  f  =  V*  "  n3ff3  r3^ 

(4  3) 

or  combining  (42)  and  (43): 


(1  -  n.  *.) 

'i  (1  -  11  -  "Vi1  =  ¥'  -  n2  Vz1 


¥' 


V  +  n 


V>  "*> 


2  (1  -  n2  rji 


-V-  X 


(1  -  n3  *3r3) 


(44) 


which  represent  more  inclusive  optimization  conditions.  -Comparison  of  (41) 
with  (44)  indicates  the  significant  difference  between  previous  '  optimisation" 
and  the  more  realistic  optimization  criteria. 


